JOURNAL OF RESEARCH of the National Bureau of Standards— D. Radio Propagation 
Vol. 66D, No. 3, May-June 1962 

Statistical Distribution of the Amplitude and Phase of 

a Multiply Scattered Field 

Petr Beckmann 

(Received December 26, 1961) 

Contribution from the Institute of Radio Engineering and Electronics, Czechoslovak Academy of Sciences, 

Prague 

The probability distribution of the amplitude and phase of the sum of a large number 
of random two-dimensional vectors is derived under the following general conditions: Both 
the amplitudes and the phases of the component vectors are random, the distributions 
being arbitrary within the validity of the Central Limit Theorem; in particular, the dis- 
tributions of tiie individual vectors need not be identical, the amplitude and phase of each 
component vector need not be independent and the distributions need not be symmetrical. 
The distributions formerly derived by Rayleigh, Rice, Hoyt, and Beckmann are shown 
to be special cases of this distribution. 



1. Introduction 

The olectrojiuiguetic field scattered by randomly 
distributed scatterers necessarily always consists 
of the individually scattered waves which mutually 
interfere to form the resulting total field. We may 
thus write 



E=re'^=^Aje'^J 



j=i 



(1) 



where r is the amplitude and 9 the phase of the 
resulting field and Aj exp (i(l>j) are the elementary 
component waves. The amplitudes Aj and phases 
(/)y, and hence also r and 6, are random quantities. 
The sum (1) may be represented in the complex 
plane as the sum of randojn vectors (fig. 1) ; the 




Figure 1. Random vector sum in the complex plane. 



problem is also identical with tlie '^random-walk 
problem'^ in inathematical statistics, for our task 
is to determine the prohal)ility distributions of r 
and 6 if tlie distributions of tJie Aj and </)j are knowii. 
Tlie problem in its most elementary fonn, wlien 
tlie Aj are constant and the 0y are all distributed 
uniformly over an interval of length 27r, was solved 
by Rayleigh [1896] and leads to the well-known 
Ray](M<r]i dislribution 



2){r)-. 



2r _,2 

-—e 
n 



T^/n 



(2) 



Other special cases of the distribution of E were 
derived by Rice [1944, 1945] and by Hoyt [1947]. 

A more general solution, namely for all Aj con- 
stant and arbitrary symmetrical distributions of 
the phases was derived by Beckmann [1959]. 
The Rayleigh, Rice, and Hoyt distributions are 
included in this distribution as special cases. 

The purpose of the present paper is to derive 
the distributions of r and 6 in the general case 
when both the Aj and the 0^ are random and 
arbitrarily distributed. The distributions of the 
individual Aj or (pj need not be identical, Aj and 
(l)j may be correlated and the distributions of tJie 
(fyj need not be symmetrical. Tlie only restriction 
is that the distributions of tlic quantities Aj cos </)^- 
and Aj sin </>,• satisfy the conditions of the Central 
Limit Theorem [GnedeidvO and Kohnogorov, 1949]. 
In physical practice these conditions are practi- 
cally always satisfied provided that n, the number 
of interfering waves, is large and that these waves 
are statistically independent.^ 

Our problem, then, is to determine the probability 
distribution p(r) and the (more rarely required) 
distribution p{d) when the two-dimensional distri- 
butions tVj(A,(l)) are given. 

1 Generalization for dependent variables is also possible [Bernstein, 1944]. 
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2. General Solution 
We introduce the quantities 

x=r cos 6] y=r sin 6 



(3) 



and denote the mean value by angular brackets 
( ) and the variance by the symbol D. Then 

n f °° r^- 
{x)=a=^^ Wj{A,cl))A cos dAdcj), (4) 

j=lj - CO J -<x, 

{y)=(3=J2\ w j{A,(t)) A sin (t>dAd(l>. (5) 

j=lj-cx, J-00 

Denoting the individual terms of the sums (4) and 
(5) by ^ 

aj= \Wj{A,(l))A COS (t) dAdcj), (6) 

13 j= I \Wj{A,ct))A sin ^ rfA/(^, (7) 

we have 

D{x}=s,^t] r I U,-(^,<^)^' cos' <l>dAdct>-af\, (8) 
^{2/1 =^2-=S [ rjw/^,«)^'sin2(^ (Mr/0-^n. (g) 

The covariance of x and y is, after elementary 
manipulations, 

cov (.r,7/)=| S ^^^,•(^,0)^' sin 2(i>dAd4> 

+ iZlZ' \ [wi{A,<t>)A cos rf^rf0 

rh/;,U,</>)^ sin <^ dAd<t>-aff, (10) 

where the apostrophe at the second sum indicates 
that the terms with i=j are to be omitted. 

We now temporarily assume that the distribution 
of is symmetrical about its mean value zero. Then 
the integrals in (5), (7), and (10) vanish owing to 
the factor sin or sin 20 in the integrands, so that 



l3j=fi=G0Y (x,y)=0. 



(11) 



Now according to the Central Limit Theorem the 
distribution of the sum of many independent random 
variables tends under certain conditions that in 
practical applications are almost invariably satisfied, 
to the normal (Gaussian) distribution.^ 



2 To save space, we henceforward omit the integration limits — oo, oo; unless 
otherwise indicated, all integrals in this paper are therefore to be taken from 

— c» to CO. 

3 For an exact statement of the Central Limit Theorem and the weakest (Linde- 
berg) conditions of its validity, see [Gnedenko and Kolmogorov, 1949, Ch. 5, 
Sec. 26], or [Richter, 1956, Ch. VII, Sec. 4]. An unnecessarily strong condition 
for the validity of the Central Limit Theorem is e.g. that the distributions of the 
terms be identical and their variance exist (Bernoulli condition). 



Since by (11) x and y are uncorrelated (and hence, 
as normal random variables, independent), the two- 
dimensional probability density of x and y is 



W{x,yy- 



27rV'§i'S2 



1 r {x-aY V'l .^^, 



Transforming to polar coordinates as in (3), we 
find the required distribution in the form 



p{r) = 



2 Try 6* 



= exp - 

1S2J0 L 



(rcos— a)^ r^si 



2<9i 



W-"]*. 



(13) 



which after elementary manipulations may be 
written as 



p{r)=Se- 



n2T 



p—Pcos^d+Q ^^^0 Jn 



(14) 



where S, T, P, and Q are functions of r, a, Su and ^2. 
The integral in (14) may be expressed as a series of 
modified Bessel functions of the first kind (1^) as 
follows [Beckmann and Schmelovsky, 1958]: 



r"^-Pcos2^+^cos^^^^_27r6-^/2^ (-l)"e^/. ^^ 

Jo Tn=0 



for m=0, 
2 for m?^0. 



where 

so that (13) may be expressed as follows: 
p{r)- 



(J)h„XQ) 

(15) 

(16) 



V-^i^: 



e 2Si iSiS2 



X 



S<-""-'"(S''')^'-("?;') <"' 



Before normalizing the distribution (15), we calcu- 
late the mean square of r, which, apart from a con- 
stant factor (l/1207r in MKS units), equals the mean 
scattered power. From (3) and the formula for the 
variance of a random quantity 

D{r}=(r}-(ry 
and from the independence of x and y it follows that 

(r'}=s,+S2+a'. (18) 

The RMS value of the scattered field equals the 
square root of this expression. 

It is easily seen from (12) that the resulting vector 

r exp (id) is the sum of a constant vector a directed 
along the x-axis and a random vector H (fig. 2); the 
X and y components of this random (Hoyt) vector 
are normally distributed with mean values zero and 
unequal variances Si and S2. 
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Figure 2. Components of the resultant vector r and its equi- 
probahility curves. 



We "now normalize in accordance with (18) bv 
using Itiie ratio of the RMS values of the constant 
and random cojuponents 



W(x^) = const 




Figure 3. Transformation of coordinates when w((/>) is not 
an even function. 



■"WV^ 
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,_L- 




B-- 



VS1+S2 

and introducing the "normalized ainijliludo" 

r 



V«l+«2 

and the "asymmetry factor" 



K- 



(1.9) 



(20) 



(21) 



Using this normahzation,* we ohtaiii tiio r(K[iiircd 
distribution by substituting (19) to (21) in (17) in 
the general form 

Xg(-l)-6./^(^^^p2)/,^[^(l+K2)p^ (22) 

The distribution (22) is determined through the 
parameters B and iiC according to (19) to (21), which 
in turn are determined through a, s^, and So accordins: 
to (4), (8), and (9). 

The probabihty densities p(p;B,K) according to 
(22) for various values of B and K are shown in 
figures 5 to 9. 

The complement of the distribution function, i.e., 



P{p>z;B,K}=l-^j^p(p;B,K)dp 



(23) 



Figure 4. Scattering by rough or turbulent layers in the 
troposphere. 



' This normalization differs from that used by Beckmann, [1959 and I960]. 




Figure 5. Probability densities p(p) for B = 0. 
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Figure 6. Probability densities p (p) for B = 0.5. 
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Figure 7. Probability densities p(p) for B=l. 
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Figure 9. Probability densities p(p) for B = 3. 



234 



for various values of B ajul K arc plotted on proba- 
bility paper in figures 10 to 14.^ 



5 The curves in figs. 5 to 14 were calculated by direct numerical integration of 
(13) by puiichod cnrd ninrliine nndor the guidan.io of Mr. K. Vicb. 
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Figure 10. Distribution curves o/ P|p>z) for K2= 1. 
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Figure 11. Distribution curves o/P{p>zt /or K2 = 2. 
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Figure 12. Distribution curves o/ Pj p>z} for K- = 3. 



As may be seen from figures 10 to 14, the distri- 
bution becomes practically^ normal for B>3, K<,Z^ 
as may also be shown theoretically (cf. appendix). 
The mean value and variance in this case is given by 



{p)^B, D{p} 



1 



From (IS), (19), and (20) we have 



<P^)=l+i^-; pkm8=Vi + 5^- 



(24) 



(25) 



Hence 
r 



I^'rms J V.PRMS J 



-j: 



zVl+^2 



v{p;B,K)<lp- (26) 
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Fig u re 13. Distrib ution curves 0/ P { p > z j for K' = 5. 
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FiGi'RE 14. Distribution curves of P|p>z| for K^— 10. 
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Figures 15 to 19 show this relation plotted on 
Rayleigh paper (on which the Rayleigh distribution 
appears as a straight line with 45° slope). 

The less often required statistical distribution of 
the phase is found froni (13) by integrating over r 
from to 00 instead of over 6 from to 27r. After 
a somewhat tedious calculation one obtains 



p{e)-- 



Ke 



'-BKl-\-K^) 



~27r{K' cos^ 6i+sin2 



[l+(?V^6^'(l-ferf(7)], 



where 



and 



(27) 



^=^^-\2(K^ 



+K' 



cos^ S+sin^ 6) 



er{G 



■^J^^ Jo 



e-''dt. 



We now return to the case when the distributions 
of the phases </)y are not symmetrical about zero. 
The general formulas (4) to (9) are still valid, but 
neither (3 in (5) nor the covariance cov (x,y) in (10) 
will vanish, so that x and y are no longer independent 
and our derivation breaks down from (11) onward. 
In this case we proceed as follows. 

We calculate the covariance cov (x,y) according to 
(10) and determine the correlation coefficient 



^^_cov {x,y) 



(28) 



We now introduce new coordinate axes x^ and y^ 
which are turned through an angle </)o with respect to 
the original axes x,y (fig. 3). The angle <^o is so 
chosen that the quantities x^=r cos 6^ and y^ = r 
sin 6^ are uncorrelated, where 

e'=d-(i>c. (29) 

For x^ and y^ to be uncorrelated (and hence, as 
normal random variables, iudependent), it is suffi- 
cient that the two-dimensional distribution W{x\y^) 
have an axis of symmetry parallel to one of the coor- 
dinate axes. Since the curves W(x,y) =const are 
concentric ellipses with center x=a, y=p, it is there- 
fore sufficient to choose the x' and y^ axes parallel to 
the axes of the ellipses. The required angle </)o then 
follows from [Bristow 1961, p. 125] 



tan 200=- 



S1—S2 



(30) 



In this new coordinate system we then proceed as 
before; the only difference is that /5' does not vanish 
now. Instead of (13) we therefore have 

(31) 
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FiGUKE 15. Distribution curves of P{r/rRMs^z! for K"=l 
plotted on Rayleigh paper. 
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Figure 16. Distribution curves of P|r/rRMs>zl for K2=:2 
plotted on Rayleigh paper. 
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Figure 17. Distribution curves of P|r/rRMs>z) for K2 = 3 
plotted on Rayletgh paper. 
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Figure 18. Distribution curves of Pjr/rRMS>zl for 1^2 = 5 
plotted on Rayleigh paper. 



Figure 19. Disiribiition curves of P(r/rHMs>z} for K2=10, 
plotted on Rdi/teinh papei . 



wlioro tJio Jueaiiijiii- of d' , a and (^^ is cvidcjit frojii 
fi^i'in'o \). TJris iiilc^t^-ral may agaiji bo ovaluatcMl as 
a.ii ijifitiilo series of Boss(4 ruiictions (cf. appoudix), 
but for practical purposes it is usually more simple 
to perl'onn the nimuMMcal intc^^ration of (31) directly. 

3. Special Cases 

We now consider some special cases of the distribu- 
tion (22). In most cases the expression (22) will 
forjually remaiji unchanged, but the values of a, Si, 
and .V2 will vary. 

l^et the elementary scattered waves be all of the 
same kind, so that the amplitude and pluise distribu- 
tions are the same for eacli wav(\ \Ve first assume 
tlie amplitudes constant and equal to unity. Let the 
probability density of the phases be synunetrical 
about zero and equal to w{(f). Then from (4), (8), 
and (9) we have 



a=n 



'\w{4)) cos (/0 
'\w{<i)) COS^0r/(/) — — 

S2=n I w{(f)) sm^(f)(Jct). 



S\=n 



(32) 



(33) 



(34) 



Table 1 gives the values of a, Su and S2 as calcu- 
l d from (32) to (34), and also of B, K, and p in 
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accordance with (19) to (21) for the normal, uniform, 
and Simpson-distributions; the symbol sine a=(sin 
a) I a is used in the table. 



Table 1. 



w{4>) 


Normal, (0)=O, 
standard deviation a 


Uniform from —a to +a 


Simpson-distributed 
from -2a to +2a 


a 


ne-V2<r2 


n sine a 


n sinc2 a 


S\ 


\<l-e--)^ 


^ (1+sinc 2a-2 sinc^ a) 


1 (H-sine2 2a+sinc^ a) 


52 


^(1-.-) 


^(1 -sine 2a) 


|(l-sine2 2a) 


B2 


,-.2 

n 


sine2a 
1— sine2 a 


sinc^ a 
1— sinc^ a 


7^2 


coth| 


1 -sine 2a 


l-sinc2 2a 


l+sinc2a-2sinc2 a 


l+sine2 2a+sine^ a 




r2 


1-2 


ri 


P^ 


nil-e-"-) 


??(1— sine2a) 


w(l-sine^ a) 



If the phases of the elementary waves are distrib- 
uted uniforjnly over the interval (— x,7r), we obtain 
for a=T from table 1 a=0, Si=S2 = nl2, or i?=0, 
K=l, p=rl-\ln; substituting these values in (22) we 
obtain, as was to be expected, the normalized 
Ravleirfi distribution 



p{p)=2pe- 



(35) 



The Ra3deigh distribution is also obtained if is 
distributed uniformly over any interval of length 
2kT{k=l,2, . . .) or if the variance of is much 
greater than t\ in which case the distribution of is 
(with some very unrealistic exceptions) arbitrary. 
We shall call a vector whose elementary components 
Aj exp (icf)) have their phases uniformly distributed 
over the interval (— 7r,7r).a ''Rayleigh vector/^ 

The distribution of the amplitude of the sum of a 
constant vector and a Rayleigh vector was found by 
Rice [1944, 1945] and detailedly anal3^zed by Norton, 
Vogler, Mansfield, and Short [1955] and Zuhrt 
[1957]. Directing the constant vector V along the 
a:-axis, we immediately have a = y, .^^1=82=71/2, hence 
by (19) to (21) B=a/^n, <l>=r/^n, K=l and on sub- 
stituting these values in (22) we find the normalized 
Rice distribution 



p(p)=2pe-^^'^f>'^U2Bp). 



(36) 



The distribution of a vector whose x and y compo- 
nents are distributed normally with mean values zero 
and unequal variances ^1 and S2 was found by Ho^^t 
[1947]. Here we have a=B=0] on substituting in 
(22) we obtain the normalized Hoyt distribution 



p(p)'- 



K^ + 1 



K 



P exp 



IrQ^) '']^^(^^0' 



(37) 



It is evident from (12) and figure 2 that the sum 
(1) may always be represented as the sum of a con- 
stant vector a and a Hoyt vector H. 

We next consider the case in which the amplitudes 
Aj of the elementary waves are random and governed 
by the (same) probability density Wa{A). If Aj and 
4)j are independent tlien it follows from the general 
formulas (4), (8), and (9) that 



a =n(A ) p;0(0) cos (j)d(j), 



(38) 



Si=^n 



cos^ (I) d(j) (39) 

n 



(A'}jw,{c^) 

S2=n(A^) ^0(0) sin^ rf0. (40) 

If, for exaniple, the phases are uniformly dis- 
tributed over an interval of length 27r, we obtain a=0, 
1 



^i = '^2=9 7i(A^}; substituting in (17) we find 



p(r)-- 



2r 



-r^ln(A^) 



~n{A^} 



(41) 



Comparing this result with (2) this will be recog- 
nized as a Rayleigh distribution, in which the num- 
ber of components n is multiplied by the mean power ^ 
of each component. 

If a Rayleigh vector consists of components with 
different (constant or random) amplitudes Aj, then 
from (4), (8), and (9) 



a=0; s,=S2=^j:{A^y 



j=i 



(42) 



Substituting this in (17) and (27), certain proper- 
ties of a Rayleigh vector with components of unequal 
but constant amplitudes postulated by Norton, 
Vogler, Mansfield, and Short [1955] are immediately 

proved as correct (1. F{r>2} = exp ( — ^^/ S ^; 



J=i 



'0/ 



2. d distributed uniforml;^ between and 27r; 3. r 
and 6 independent). 

If the random amplitudes and phases of the ele- 
mentary waves are correlated, formulas (4), (8), 
and (9) ma^^ not be simplified; however, it is still 
true that the mean power of the random component 
equals the sum of the mean powers of the individual 
(identically distributed) elementary components, for 
in this case we have from (8) and (9) 



-§1+^2 = '^^ 



,JA'^jw{A,<t>)dcl>yiA-^^ 

=nfA'wAA)dA-^=n{A')-^^ (43) ,' 



6 More precisely the mean square of the amplitude. This differs from the mean 
power by a constant factor F, which in MKS units equals 1207r for propagation in 
free space. The distinction is immaterial for our present purposes and will be 
disregarded; if the reader objects to this procedure, he may consider (1) to present 
n sinusoidal voltages interfering across a one-ohm resistor, in which case F equals 
unity. 
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so tliat by (18) 



(r')=n{A^)+a'(l-iy 



(44) 



Since n is by assumption large, the second term 
will practically equal a^, i.e., the power of the con- 
stant component. Thus the power of the random 
component equals n{A^), or the sum of the mean 
powers of the individual components. 

It is instructive to observe the transition from a 
purely coherent field (mean power equal to n^) to a 
purely incoherent jfield (mean power equal to n). 
This depends on the phase distribution w{(l)) ; if the 
phases are constant, i.e., D{(t)}=0, then (r^) = n^, 
whereas for phase distributions with large variances, 
i.e., for !>{(/)} >>7r^ the mean power {r^)=n. Thus 
for example from (18) and table 1 we find for nor- 
mally distributed phases 



(r2)=n'6-'^'+n(l-6-'^') 



(45) 



yielding (r^)=7i^ for g=0, but (r^}=n for c-^ co (o-> >7r) . 
Similarly, for uniformly distributed phases we find 



(r'^)=n^ sinc^ a+n(l— sinc^ a) 



(46) 



again yielding {r^}=n^ for a=0, but r^=n for a-^co 
(a>>7r), and also for a=2kT. 

4. Conclusion 

The statistical distribution of the amplitude and 
phase of a multiply scattered elect roinagnetic field 
is equal to the statistical distribution of the sum of 
two-dimensional vectors with random amplitudes 
and phases. Wlien these phases are distributed 
symmetrically, the amplitude distribution of the 
resulting vector is given by (17) or in the normalized 
form (22), or by the curves of figures 5 to 19; the 
phase distribution is given by (27). In the general 
case, which includes asynnnetrical pliase distribu- 
tions, tlie resulting distribution is given by the in- 
tegral (31). Various distribution laws of the am- 
plitudes and phases of the elementary vectors change 
the values of a, Si, and S2, but not the general form 
of the above formulas. The distributions derived 
by Kayleigli [1896], Rice [1944-5], Hoyt [1947], and 
Beckmann [1959] are special cases of the above 
distribution. 

The distribution derived here is met, among other 
cases, in the propagation of radio waves in irregular 
terrain and in tropospheric scatter propagation, 
since in both cases scattering from rough surfaces is 
involved. From the above derivation it is seen that 
the amplitude of a field consisting of very many 
elementary scattered waves is not necessarily Ray- 
leigh-distributed (as is often erroneously assumed), 
but that the Rayleigh distribution, even in its most 
general form, is met only if the phases of the in- 
dividual scattered waves are distributed uniformly 
over an interval of length 27r or in some equivalent 
way indicated after equation (35). In practice this 



will not be the case if, for example, the scatterers 
are distributed in space in such a way that the 
variance of path lengths between source and point 
of observation is smaller than one wavelength. Such 
a case is shown in figure 4, where a rough or tm^bulent 
layer is assumed to be normally distributed about a 
mean level (M) with variance a and the condition 



2(x sin 7<CX 



(47) 



holds. This condition is very often satisfied in 
practice, especially for the longer wavelengths X; 
experimental measurejnents of tropospheric propaga- 
tion beyond the horizon in tJie meter band [Beck- 
mann, 1960] have in fact shown distributions as in 
figures 15 to 19 more often than a pure Rayleigh 
distribution. 

5. Appendix 

To evaluate the integral (31), one may use a result 
derived by Chytil [1961], which after elementary 
modifications reads 



r2ir 

u(P,Q,R)= exp (-P2 cos2 O+Qco^ d+R sin d)^ 

= 2Te-¥j: {-iy'eJrn(~)l2rnUWTR') 
m=0 \^ / 



Id 



COS 2m ( arctan jr j (48) 



reducing to tlu> forjnula (15) derived by Beckmann 
and Schmelovsky [1958] for E=0. For 7^=0 and 
large Q((^>P>>1) one obtains by saddle-point 
integration [Beckmann and Schmelovsky, 1958] 



where 



^(P,«,0)«/^V-S^„,(|)'" (49) 

_ 1.3.5...(2m-l) . _ , 

^'""- 2.4.6... 2m ' "^»-^- ^''^' 



Using (49) to evaluate (14), wc find after normal- 
izing by (19) to (21) for Q>2P»1. 



Pip)- 



,exp^-^-^^^(B-py'] 



K 



■yj2wB 



Now if 






_ CO / JZ2 1 \ Wl 

(52) 



this expression will obviously be negligibly small for 
all values of p except in the neighborhood of p=B, 
where the exponential factor will dominate, the 
terms with p(i/2)-m being either negligible or practically 
constant in this short interval. But the exponential 
is that of tlie normal distribution with mean value 
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{p)=B and variance Z>{p} =1/(1 +J^^); hence for 
large B (52), the distribution of p becomes normal. 
As may be seen from figures 10 to 14, in practice 
this is the case for 5>3, iP<5. 
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